Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  844 1708 (N35E24 L=196)   660    104        1 ( 806/1666- 893/1738)  56.3 E26-E20 N31-N37
      2  940 1593 (N28E14 L=206)   323     77        0 ( 917/1553- 970/1618)  71.6 E15-E12 N24-N30
      3  747 1209 (N01E19 L=202)   825    287        0 ( 687/1080- 789/1330) -53.1 E24-E14 S05-N08
      4 1603 1215 (N14W34 L=254)   293     41        1 (1576/1191-1630/1237) -33.7 W31-W35 N12-N15
      5 1902 1122 (N15W62 L=283)  2586    351        1 (1853/1034-1940/1198)  65.6 W56-W67 N08-N19
      6  560  977 (S14E28 L=192)  1082    333        0 ( 447/ 867- 688/1088)  43.0 E35-E20 S21-S06
      7  369  991 (S15E42 L=178)   832    281        4 ( 267/ 931- 473/1040)  20.6 E51-E34 S18-S11
      8 1887  967 (N05W63 L=283)   444     71        1 (1864/ 938-1907/ 989) -21.8 W59-W65 N03-N06
      9 1203  817 (S16W14 L=234)   402    168        0 (1135/ 759-1239/ 861) -35.5 W09-W16 S18-S13
     10  494  650 (S35E33 L=187)   316     65        1 ( 463/ 629- 517/ 679) -45.0 E35-E30 S36-S33
     11  861  214 (S65W02 L=222)   514     76        1 ( 835/ 190- 889/ 242)  -8.1 E01-W06 S67-S62
  

    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