Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1534 1591 (N21W53 L=209)  1284    190        3 (1455/1540-1628/1630) -14.0 W44-W61 N17-N22
      2 1145 1574 (N33W26 L=182)   665     89        1 (1099/1546-1194/1598) -26.6 W22-W29 N31-N35
      3  317 1517 (N56E54 L=102)   467     48        1 ( 291/1501- 347/1530)   0.0 E60-E47 N54-N58
      4 1758 1355 (N02W58 L=214)  1006    161        1 (1719/1270-1788/1425)  58.0 W55-W61 S03-N05
      5  691 1394 (N36E10 L=145)   221     52        0 ( 658/1380- 720/1406) -26.6 E12-E08 N35-N37
      6  150  896 (N18E70 L= 86)   946    104        3 ( 121/ 835- 178/ 934) -40.6 E75-E66 N14-N20
      7 1317  804 (S14W11 L=166)   717    263        5 (1234/ 706-1379/ 905)  36.0 W03-W15 S18-S09
      8  846  696 (S08E19 L=137)   847    262        2 ( 771/ 590- 919/ 803) -90.0 E21-E17 S15-N00
      9 1327  563 (S29W05 L=160)   337     94        0 (1274/ 547-1366/ 577) -26.6 W01-W07 S31-S27
     10  983  373 (S31E23 L=132)   882    181        1 ( 927/ 299-1022/ 434)  76.0 E26-E20 S37-S25
  

    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