Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  773 1679 (N46E25 L=127)  2883    510        8 ( 642/1603- 933/1731)   0.0 E36-E09 N39-N51
      2 1136 1629 (N44W08 L=160)   510    118        2 (1104/1583-1180/1661)  45.0 W04-W11 N40-N46
      3 1542 1537 (N37W40 L=193)   415     92        1 (1503/1517-1578/1558)  15.9 W36-W43 N35-N38
      4  770 1330 (N21E18 L=135)   504     98        2 ( 733/1290- 806/1384) -45.0 E20-E15 N18-N24
      5  320 1210 (N12E51 L=102)   613    132        0 ( 294/1154- 350/1263) -45.0 E53-E47 N08-N14
      6  708 1227 (N15E21 L=132)   878    266        5 ( 627/1184- 820/1250)   4.4 E26-E13 N11-N16
      7  810 1179 (N12E14 L=139)   286     62        1 ( 791/1149- 827/1213)  76.0 E15-E12 N10-N14
      8  667 1137 (N09E23 L=130)   681    189        1 ( 610/1095- 757/1200) -24.4 E27-E16 N06-N12
      9 1038  904 (S04W01 L=154)   326    120        0 (1022/ 848-1055/ 965)  90.0 W00-W02 S07-N00
     10  204  794 (S15E63 L= 90)  1100    211        1 ( 158/ 761- 229/ 840) -11.3 E69-E59 S17-S11
     11  393  452 (S38E52 L=101)  1490    374        1 ( 338/ 399- 449/ 507)  30.5 E62-E44 S43-S33
     12 1140  381 (S40W11 L=164)  1206    292        1 (1007/ 343-1313/ 426)  14.0 W00-W24 S43-S36
  

    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