Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1471 1736 (N58W38 L=171)  1172    163        2 (1395/1717-1549/1753)  11.3 W25-W56 N54-N60
      2  917 1641 (N40E18 L=115)   324     83        2 ( 877/1606- 944/1662) -26.6 E21-E15 N37-N41
      3 1426 1624 (N46W26 L=159)   314     44        0 (1397/1608-1452/1636)  18.4 W23-W29 N45-N47
      4  655 1452 (N23E31 L=102)   488    152        1 ( 615/1396- 676/1519) -77.5 E33-E29 N18-N27
      5 1205 1463 (N31W07 L=139)   246     85        2 (1187/1425-1220/1509) -68.2 W05-W07 N28-N33
      6  480 1423 (N19E44 L= 89)   680    144        3 ( 437/1393- 542/1468) -16.7 E48-E38 N16-N21
      7  575 1407 (N19E36 L= 97)   349     90        1 ( 559/1364- 593/1446) -68.2 E37-E33 N16-N21
      8 1738 1415 (N34W55 L=187)   948    154        1 (1671/1371-1792/1446) -13.0 W47-W61 N31-N35
      9  455  817 (S19E36 L= 97)   845    151        2 ( 413/ 754- 509/ 857)  41.2 E39-E31 S23-S15
     10 1702  650 (S14W55 L=187)   396     30        1 (1684/ 629-1725/ 669)  33.7 W52-W56 S15-S13
     11 1400  411 (S34W39 L=172)   648    193        1 (1328/ 363-1488/ 453)  41.6 W34-W44 S38-S30
  

    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