Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  811 1783 (N35E43 L=308)   330     39        1 ( 792/1761- 832/1803) -18.4 E45-E39 N33-N36
      2 1098 1697 (N37E18 L=332)   299     43        0 (1073/1673-1126/1717)  -9.5 E20-E14 N36-N38
      3  606 1556 (N14E42 L=308)   275     63        0 ( 580/1538- 635/1577)   0.0 E44-E39 N13-N15
      4  298 1173 (S15E52 L=298)  1092    160        1 ( 268/1094- 334/1253)  84.3 E54-E49 S19-S09
      5 1700 1189 (N24W38 L= 28)   337     79        0 (1662/1166-1741/1204)   7.1 W34-W42 N22-N25
      6 1006  805 (S18W05 L=355)   203     73        0 ( 988/ 768-1029/ 835)  90.0 W04-W05 S20-S15
      7 1257  731 (S15W22 L= 12)   265     46        0 (1227/ 707-1290/ 753)  63.4 W20-W23 S16-S12
      8 1334  617 (S19W31 L= 21)   471    161        0 (1257/ 598-1415/ 634)  33.7 W26-W35 S21-S15
      9  776  579 (S38E02 L=348)   760    162        2 ( 701/ 529- 837/ 616) -10.3 E08-W03 S39-S35
     10 1410  370 (S30W51 L= 41)   827    106        1 (1367/ 332-1462/ 402)  45.0 W46-W53 S32-S26
  

    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