Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  972 1435 (N27W07 L=335)   351    115        0 ( 919/1402-1023/1470) -56.3 W03-W08 N24-N30
      2 1817 1274 (S02W59 L= 28)  3562    481        6 (1664/1077-1879/1368) -59.5 W46-W69 S13-N04
      3 1334 1296 (N11W24 L=352)   322     96        0 (1285/1274-1378/1314) -38.7 W21-W26 N09-N13
      4 1504  963 (S11W27 L=355)   310     59        0 (1489/ 918-1521/ 993)  63.4 W25-W28 S13-S09
      5 1342  697 (S23W12 L=340)  2710    787        2 (1002/ 523-1631/ 863)   7.6 E12-W33 S26-S19
      6 1072  799 (S11E02 L=326)   463    166        0 ( 990/ 771-1140/ 845)   0.0 E07-W02 S11-S09
      7  210  798 (N06E61 L=267)   791    117        0 ( 171/ 763- 266/ 834) -32.0 E66-E56 N03-N08
      8  149  763 (N05E72 L=256)   765     65        0 ( 125/ 729- 176/ 792)  15.3 E78-E66 N02-N06
      9  591  536 (S16E39 L=290)   281     45        2 ( 560/ 519- 612/ 558) -33.7 E41-E36 S17-S14
     10  952  447 (S30E19 L=309)   350    107        0 ( 931/ 394- 985/ 485)  29.7 E22-E15 S32-S27
     11  495  285 (S30E66 L=262)   645     41        1 ( 473/ 260- 521/ 310)   8.7 E72-E59 S32-S28
  

    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