Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1406 1634 (N44W33 L=131)  1564    298        7 (1346/1565-1455/1691)  29.7 W26-W40 N38-N48
      2  934 1477 (N32E07 L= 91)   324     90        1 ( 920/1437- 957/1522) -80.5 E08-E05 N29-N35
      3 1727 1462 (N30W56 L=154)   481     85        1 (1701/1439-1751/1480) -26.6 W54-W58 N28-N31
      4  575 1320 (N21E31 L= 66)   470    109        2 ( 554/1261- 593/1369) -81.9 E32-E29 N17-N24
      5 1556 1312 (N21W37 L=134)   514    117        2 (1510/1269-1611/1348) -29.7 W33-W40 N18-N23
      6 1556 1226 (N15W36 L=133)   582    175        2 (1488/1202-1620/1260)   5.2 W30-W41 N14-N17
      7 1103 1164 (N12W05 L=102)   366     95        0 (1081/1118-1132/1219) -63.4 W03-W06 N09-N15
      8  595 1131 (N10E28 L= 70)   225     79        0 ( 562/1111- 634/1153) -18.4 E30-E24 N08-N10
      9 1565 1096 (N07W36 L=133)   394    109        0 (1508/1073-1613/1128) -20.6 W31-W39 N05-N09
     10  494  991 (N01E35 L= 63)   483    137        0 ( 444/ 949- 553/1030)  26.6 E38-E30 S02-N03
     11  393  697 (S18E45 L= 53)   343     86        1 ( 380/ 651- 406/ 732)  71.6 E45-E43 S21-S15
     12  721  625 (S22E20 L= 77)   653    138        1 ( 665/ 600- 791/ 659) -14.0 E23-E15 S23-S19
     13 1578  614 (S23W40 L=137)   510     75        2 (1543/ 585-1614/ 639)  38.7 W36-W42 S25-S21
     14  487  592 (S25E39 L= 59)   574    158        1 ( 436/ 573- 551/ 612)  12.5 E42-E33 S26-S23
     15 1117  405 (S38W08 L=105)   561    107        2 (1087/ 363-1146/ 449)  90.0 W05-W09 S41-S34
     16  989  347 (S43E03 L= 95)   974    156        1 ( 934/ 315-1035/ 412)  58.0 E07-W01 S45-S37
     17  720  375 (S41E25 L= 72)   404     76        0 ( 694/ 353- 740/ 408)  45.0 E27-E23 S42-S38
  

    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