Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  508 1494 (N31E42 L= 37)   713    118        2 ( 449/1470- 568/1513)  -4.8 E48-E36 N29-N32
      2  925 1362 (N24E08 L= 71)   453    152        0 ( 855/1350-1006/1378)  -5.7 E12-E02 N23-N24
      3 1493 1201 (N15W31 L=109)   297    114        0 (1438/1181-1532/1217)  -8.1 W26-W33 N13-N15
      4  162  861 (S11E69 L=  9)  2153    228        6 ( 122/ 756- 191/ 955) -81.3 E76-E65 S18-S05
      5 1417  839 (S08W26 L=104)   266     54        0 (1391/ 811-1441/ 863)  45.0 W23-W27 S09-S06
      6 1659  735 (S14W46 L=124)  1197    206        3 (1596/ 671-1725/ 819)  45.0 W41-W50 S18-S08
      7  378  698 (S20E46 L= 32)  1008    190        1 ( 346/ 614- 430/ 782) -74.7 E48-E42 S25-S14
      8 1070  493 (S31W05 L= 84)   559    126        2 (1051/ 428-1097/ 552) -66.0 W03-W07 S36-S27
      9 1118  375 (S41W10 L= 89)   365     59        0 (1093/ 347-1144/ 405)  53.1 W07-W11 S43-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