Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  648 1395 (N33E15 L=178)  1156    242        2 ( 544/1325- 726/1439) -31.6 E23-E10 N27-N37
      2 1540 1322 (N08W38 L=232)   552    139        3 (1508/1284-1594/1351) -31.0 W35-W41 N05-N09
      3  148 1064 (N23E62 L=132)  1215    198        2 ( 112/ 982- 173/1121) -66.0 E65-E59 N18-N27
      4  840  783 (S06E16 L=178)  3395    503        9 ( 663/ 566- 985/ 936)  -8.4 E26-E08 S20-N03
      5 1247  863 (S11W09 L=203)   374    201        0 (1179/ 845-1328/ 878) -20.6 W05-W13 S13-S09
      6 1265  757 (S17W08 L=202)   297    142        1 (1206/ 723-1341/ 787)   0.0 W03-W12 S18-S16
      7  667  645 (S11E30 L=164)   439    113        1 ( 633/ 598- 703/ 688)  29.7 E33-E26 S13-S08
      8 1049  555 (S25E10 L=184)   525    140        0 ( 981/ 520-1130/ 598)   0.0 E14-E03 S25-S23
      9  803  525 (S21E26 L=168)   492    112        2 ( 754/ 499- 854/ 555)   0.0 E29-E21 S21-S18
     10  511  300 (S29E58 L=135)   501     89        1 ( 484/ 267- 532/ 339)  71.6 E60-E56 S31-S25
  

    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