Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1670 1481 (N27W54 L= 80)   839    180        1 (1611/1451-1722/1503)  -5.2 W48-W59 N24-N28
      2 1372 1299 (N19W25 L= 51)   499    106        3 (1318/1265-1424/1331) -23.2 W21-W28 N16-N21
      3  278 1168 (N17E53 L=333)   679    163        3 ( 238/1140- 320/1202) -20.6 E57-E49 N15-N19
      4  646 1147 (N15E24 L=  3)   317    102        0 ( 632/1091- 660/1191) -71.6 E24-E22 N11-N17
      5 1052 1075 (N08W02 L= 28)   228     53        1 (1027/1043-1079/1106) -76.0 W00-W03 N05-N09
      6  172 1011 (N07E66 L=320)  2095    316        4 ( 136/ 940- 208/1079) -76.0 E72-E61 N03-N11
      7 1787  905 (S10W54 L= 80)  1024    140        3 (1765/ 841-1806/ 982) -83.7 W51-W55 S14-S05
      8 1828  781 (S19W59 L= 85)  1966    328        4 (1771/ 734-1902/ 848)   7.1 W52-W68 S21-S14
      9 1406  671 (S20W23 L= 49)   256     80        0 (1368/ 640-1433/ 716)  38.7 W19-W24 S22-S17
     10  456  537 (S23E45 L=341)   571    111        2 ( 444/ 489- 474/ 590)  60.3 E48-E43 S26-S19
  

    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