Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  411 1401 (N18E49 L=242)   945    220        2 ( 380/1303- 437/1493) -53.1 E53-E43 N12-N24
      2  879 1317 (N18E12 L=278)   272    150        1 ( 836/1274- 922/1362) -26.6 E15-E09 N15-N20
      3  815 1048 (N01E13 L=277)   244     67        1 ( 786/1022- 849/1071)  33.7 E14-E11 N00-N02
      4  839  859 (S10E10 L=280)   476    110        2 ( 812/ 818- 867/ 890)  53.1 E11-E08 S12-S08
      5  587  798 (S17E26 L=264)   246     48        0 ( 565/ 777- 605/ 831)  90.0 E28-E25 S18-S14
      6  875  753 (S17E07 L=284)   535    185        1 ( 791/ 732- 941/ 798)  -5.7 E12-E02 S17-S14
      7 1538  716 (S13W38 L=328)   840    191        0 (1509/ 642-1585/ 779)  83.7 W35-W40 S17-S08
      8 1303  700 (S16W21 L=312)   492    194        1 (1234/ 657-1354/ 772)  58.0 W17-W23 S19-S11
      9  211  741 (S25E65 L=226)   596     84        1 ( 173/ 729- 241/ 769)  15.9 E73-E59 S26-S22
     10  727  547 (S32E16 L=274)   324    105        1 ( 676/ 533- 772/ 568)   0.0 E20-E12 S33-S30
     11 1483  359 (S38W49 L=340)  1068    105        2 (1440/ 324-1522/ 407)  80.5 W44-W53 S41-S34
  

    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