Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1460 1313 (N13W33 L= 57)   231     89        0 (1426/1281-1496/1345) -71.6 W32-W35 N09-N15
      2  491 1264 (N32E31 L=353)   268     93        0 ( 446/1232- 533/1286) -29.7 E34-E27 N29-N34
      3 1119 1235 (N17W11 L= 34)   618    189        0 (1023/1203-1198/1271) -37.9 W05-W14 N13-N20
      4  835 1203 (N22E07 L= 16)   356    166        0 ( 763/1180- 916/1233) -33.7 E11-E02 N18-N24
      5  687  946 (N10E21 L=  2)   539    244        2 ( 632/ 841- 727/1050)  85.6 E23-E19 N03-N16
      6 1673  881 (S17W37 L= 61)   498    130        1 (1638/ 858-1734/ 919) -35.5 W34-W42 S19-S14
      7 1129  709 (S13E01 L= 22)   199     64        1 (1098/ 685-1158/ 734) -76.0 E02-W00 S15-S11
      8 1734  571 (S41W44 L= 67)  1791    302        2 (1627/ 479-1849/ 647)   8.1 W31-W60 S44-S37
      9  604  514 (S14E39 L=344)   985    147        2 ( 571/ 433- 648/ 571)  52.1 E42-E34 S19-S10
     10 1448  384 (S45W14 L= 37)   476     81        1 (1401/ 361-1491/ 408) -14.0 W09-W17 S47-S43
     11  522  264 (S31E67 L=317)   807     58        1 ( 504/ 234- 546/ 293)  15.5 E78-E60 S33-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