Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  671 1579 (N35E23 L=325)   276     50        1 ( 648/1557- 693/1604)  56.3 E24-E20 N33-N36
      2  650 1412 (N23E23 L=325)   245     64        0 ( 633/1382- 673/1448) -63.4 E23-E20 N21-N25
      3 1697 1351 (N15W48 L= 35)   549    120        2 (1652/1310-1746/1407) -45.0 W44-W51 N12-N18
      4 1474 1272 (N11W30 L= 17)   332     89        1 (1453/1239-1503/1322) -59.0 W28-W31 N08-N13
      5 1044 1172 (N06W02 L=349)   973    260        2 ( 934/1062-1120/1255) -50.7 E04-W06 N00-N11
      6  653 1056 (N02E22 L=325)   403    162        1 ( 597/1032- 738/1093) -26.6 E25-E17 N00-N04
      7 1659  843 (S16W41 L= 29)   536    104        3 (1626/ 781-1688/ 897)  66.8 W39-W43 S19-S12
      8  194  861 (S07E61 L=287)   657    150        2 ( 161/ 820- 234/ 891)  18.4 E66-E56 S08-S05
      9  862  762 (S17E11 L=336)   237     89        1 ( 823/ 733- 890/ 791) -45.0 E13-E09 S19-S15
     10 1699  580 (S32W49 L= 36)   353     67        2 (1672/ 564-1730/ 595)   0.0 W46-W51 S33-S31
     11 1434  446 (S41W29 L= 17)  1446    346        3 (1282/ 372-1589/ 521)  16.3 W17-W41 S46-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