Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1718 1357 (N19W54 L= 33)   618    117        2 (1668/1329-1765/1385) -26.6 W49-W58 N17-N21
      2  220 1243 (N20E62 L=277)   746    117        2 ( 196/1194- 252/1296) -36.9 E67-E57 N17-N23
      3  418 1067 (N09E40 L=299)   274    110        0 ( 405/1009- 443/1106)  71.6 E41-E38 N05-N11
      4  576  991 (N04E29 L=310)   813    221        2 ( 543/ 889- 604/1093)  74.7 E31-E26 S01-N10
      5  273  895 (S01E54 L=285)  3208    500        3 ( 218/ 727- 415/1080) -51.3 E59-E42 S11-N10
      6 1024  865 (S05E01 L=338)   343     85        0 ( 981/ 843-1061/ 896)  26.6 E03-W01 S06-S03
      7 1245  733 (S15W12 L=351)   433    107        1 (1192/ 699-1285/ 760) -38.7 W09-W14 S17-S12
      8  469  686 (S14E40 L=299)   334     52        1 ( 434/ 664- 496/ 702)   0.0 E42-E37 S15-S13
      9  464  580 (S21E43 L=296)   829    188        0 ( 381/ 542- 546/ 625) -26.6 E49-E36 S24-S18
     10  657  344 (S40E36 L=303)  1680    280        0 ( 569/ 245- 745/ 448)  26.6 E58-E24 S49-S31
  

    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