Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  770 1734 (N59E11 L= 50)  1133    279        3 ( 645/1695- 898/1770)   0.0 E26-W04 N55-N60
      2 1044 1436 (N31W08 L= 68)   332    103        1 ( 994/1394-1074/1459) -51.3 W04-W09 N28-N33
      3  825 1368 (N30E08 L= 52)   315    106        0 ( 787/1341- 876/1405) -45.0 E10-E05 N27-N32
      4  236 1198 (N25E59 L=  2)   662    167        2 ( 201/1141- 281/1271) -33.7 E65-E52 N21-N29
      5 1546 1236 (N11W37 L= 97)   337     70        1 (1512/1212-1577/1261) -36.9 W34-W38 N08-N12
      6  394 1128 (N19E42 L= 19)   365    112        1 ( 381/1080- 416/1173) -63.4 E42-E39 N16-N22
      7 1908  965 (S13W68 L=128)  2267    210        4 (1879/ 891-1941/1060) -74.7 W63-W75 S18-S07
      8 1542  954 (S06W31 L= 92)   261     69        0 (1530/ 915-1554/ 995) -90.0 W30-W32 S08-S03
      9  169  957 (N10E66 L=354)   764     56        1 ( 141/ 924- 195/ 981)  31.0 E70-E62 N08-N11
     10 1200  744 (S13W07 L= 67)   352    110        0 (1158/ 719-1253/ 777)  26.6 W04-W10 S14-S11
     11 1764  630 (S32W52 L=113)  1484    203        0 (1682/ 561-1851/ 683)   8.1 W43-W65 S35-S29
     12  381  589 (S14E52 L=  8)   714    178        1 ( 343/ 508- 401/ 663)  42.5 E60-E48 S19-S08
     13  239  604 (S12E69 L=352)   646     56        1 ( 212/ 577- 260/ 623) -26.6 E73-E65 S14-S11
     14 1269  360 (S42W08 L= 68)   505     94        1 (1228/ 333-1322/ 379)   0.0 W04-W12 S43-S39
     15 1080  266 (S47E10 L= 51)   994    129        1 (1031/ 230-1148/ 291)   9.5 E15-E03 S49-S44
     16  908  202 (S50E29 L= 31)  1241    137        1 ( 869/ 158- 947/ 257)  33.7 E38-E22 S55-S44
  

    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