Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  534 1646 (N43E48 L= 42)   364     56        0 ( 516/1627- 555/1670) -20.6 E52-E44 N41-N44
      2 1416 1570 (N40W31 L=122)   357     41        1 (1391/1544-1443/1594) -36.9 W29-W33 N37-N41
      3  300 1516 (N31E69 L= 22)   983    135        3 ( 279/1448- 327/1577) -22.8 E79-E59 N27-N35
      4 1681 1452 (N31W53 L=144)   624    110        1 (1639/1409-1717/1487) -38.7 W49-W56 N27-N33
      5 1219 1366 (N25W13 L=103)   410    127        1 (1170/1326-1271/1403)  35.5 W09-W16 N22-N27
      6 1496 1238 (N17W31 L=122)   999    325        2 (1390/1136-1614/1328) -38.7 W24-W39 N10-N22
      7 1855 1226 (N16W67 L=158)  2065    254        3 (1827/1129-1898/1313) -90.0 W61-W72 N09-N21
      8  366 1138 (N08E46 L= 45)   424    200        0 ( 308/1057- 393/1186) -48.8 E51-E43 N02-N10
      9 1750  855 (S07W52 L=143)   528     41        1 (1729/ 828-1775/ 880)  71.6 W50-W54 S08-S05
     10  905  738 (S15E07 L= 84)  1267    317        4 ( 849/ 617- 952/ 875)  83.3 E10-E03 S23-S06
     11 1058  836 (S08W02 L= 93)   219    109        0 (1016/ 820-1092/ 864)  26.6 W00-W04 S09-S06
     12  459  785 (S14E38 L= 53)   305     52        0 ( 434/ 771- 493/ 806)  11.3 E40-E35 S14-S12
     13  186  743 (S19E70 L= 21)   655     51        0 ( 158/ 720- 218/ 760)  11.3 E79-E63 S20-S17
     14  472  733 (S17E37 L= 53)   276     67        2 ( 450/ 708- 496/ 755)  45.0 E39-E35 S18-S15
     15 1207  672 (S19W13 L=104)   347    113        1 (1172/ 625-1238/ 712) -45.0 W10-W15 S21-S16
     16 1476  528 (S28W35 L=126)  1634    240        2 (1436/ 441-1520/ 652) -63.4 W31-W42 S35-S19
     17  801  597 (S25E14 L= 77)   392    113        1 ( 734/ 581- 842/ 623)   7.1 E19-E11 S25-S22
     18 1680  495 (S31W58 L=149)   587     93        1 (1637/ 483-1723/ 511)  -8.7 W51-W64 S32-S29
  

    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