Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  492 1682 (N63E31 L= 39)   488     48        0 ( 466/1659- 516/1699) -10.3 E37-E26 N60-N65
      2  877 1521 (N37W05 L= 74)   445    108        0 ( 820/1495- 933/1541) -45.0 W01-W07 N33-N39
      3  483 1483 (N46E28 L= 42)   849    154        1 ( 408/1455- 561/1504)  -7.1 E36-E20 N43-N47
      4  444 1197 (N28E33 L= 36)   495    143        1 ( 420/1149- 483/1260) -66.8 E34-E30 N24-N31
      5 1537 1020 (S09W30 L= 99)   242     78        0 (1517/ 987-1554/1060)  53.1 W27-W31 S10-S06
      6  639 1025 (N13E22 L= 47)   332    100        0 ( 611/ 985- 679/1052)  21.8 E24-E19 N11-N14
      7 1513  752 (S24W23 L= 92)   712    236        2 (1480/ 678-1566/ 872) -90.0 W20-W25 S29-S16
      8 1705  695 (S34W38 L=108)   517     72        1 (1685/ 649-1722/ 739) -74.1 W36-W40 S37-S30
      9 1314  640 (S25W07 L= 77)   270     56        0 (1277/ 628-1350/ 651) -21.8 W04-W09 S26-S23
     10  625  316 (S28E51 L= 19)  2504    313        1 ( 569/ 196- 743/ 426)  18.4 E73-E34 S36-S23
     11 1460  346 (S53W13 L= 82)   426     54        0 (1431/ 326-1490/ 365) -39.8 W10-W16 S55-S49
     12 1319  236 (S58E05 L= 65)   623     58        0 (1297/ 214-1349/ 263) -66.8 E07-E01 S60-S53
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y