Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  869 1724 (N31E34 L= 30)   339     72        0 ( 839/1693- 892/1754) -80.5 E35-E32 N27-N33
      2  563 1678 (N20E53 L= 11)   355     56        0 ( 533/1656- 587/1697)  56.3 E54-E50 N17-N21
      3  372 1492 (N04E58 L=  7)   802    124        2 ( 322/1449- 406/1554)  76.0 E60-E54 N00-N08
      4  735 1332 (N03E25 L= 39)   384    159        1 ( 703/1273- 776/1405) -36.9 E28-E20 N00-N06
      5  986 1096 (S04E04 L= 60)   504    175        0 ( 935/1027-1029/1172)  90.0 E05-E02 S08-N01
      6  255  960 (S29E53 L= 11)   650    105        1 ( 219/ 926- 280/1010)  45.0 E57-E49 S32-S26
      7  996  827 (S18W04 L= 68)   661    217        3 ( 910/ 776-1087/ 886)  59.0 W00-W08 S22-S12
      8  759  847 (S23E10 L= 54)   276    169        0 ( 690/ 804- 832/ 867)  11.3 E14-E04 S25-S22
      9  497  784 (S34E28 L= 37)   250     74        0 ( 463/ 769- 537/ 801)  26.6 E30-E24 S35-S32
     10  591  650 (S40E17 L= 47)   992    248        0 ( 502/ 583- 698/ 727)  49.1 E24-E11 S47-S32
     11 1564  378 (S24W64 L=128)  1529    145        3 (1501/ 353-1625/ 424)  37.9 W55-W72 S27-S20
     12 1226  216 (S46W56 L=120)  1661    171        3 (1190/ 162-1285/ 295)  -8.7 W44-W70 S47-S41
  

    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