Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  231 1369 (N28E68 L=149)   538     59        1 ( 204/1344- 264/1387)  14.0 E72-E62 N26-N29
      2  826 1169 (N15E12 L=205)   283    102        0 ( 778/1156- 872/1183)   0.0 E14-E08 N13-N15
      3 1168  894 (S04W08 L=225)   597    231        2 (1072/ 860-1271/ 938) -22.6 W02-W14 S06-N00
      4  188  776 (S09E68 L=148)   869     85        2 ( 166/ 742- 210/ 819) -45.0 E71-E63 S11-S06
      5 1613  674 (S22W40 L=256)   541    162        1 (1562/ 619-1656/ 737)  49.4 W36-W42 S25-S18
      6 1548  658 (S22W35 L=251)   298     92        0 (1524/ 620-1570/ 707)  59.0 W32-W35 S24-S19
      7  337  633 (S18E54 L=162)   850    163        2 ( 291/ 580- 382/ 706) -76.0 E57-E50 S21-S13
      8  609  610 (S20E31 L=186)   783    175        5 ( 561/ 546- 677/ 674)  28.6 E36-E25 S23-S15
      9  530  446 (S31E43 L=173)  1397    238        3 ( 512/ 319- 561/ 544)  54.2 E53-E37 S42-S24
     10 1600  482 (S36W44 L=261)   587     59        2 (1576/ 449-1626/ 510)  76.0 W41-W46 S38-S33
  

    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