Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1616 1577 (N18W58 L=223)  3128    625        7 (1475/1539-1779/1608) -22.8 W45-W78 N11-N24
      2 1500 1363 (N11W38 L=204)   296    108        0 (1486/1303-1515/1411)  56.3 W36-W40 N07-N13
      3 1897 1205 (S11W67 L=233)  1235    183        2 (1873/1139-1917/1268)  54.5 W62-W71 S14-S07
      4  592 1174 (N26E22 L=144)   820    255        1 ( 498/1074- 665/1246) -63.4 E28-E19 N18-N32
      5 1837  961 (S22W52 L=218)  3087    359        2 (1782/ 795-1866/1125)  84.6 W46-W55 S31-S10
      6  139  929 (N21E67 L= 99)   739    132        0 ( 129/ 869- 147/ 996)  58.0 E69-E64 N16-N24
      7 1691  786 (S27W35 L=201)  1189    267        1 (1659/ 681-1736/ 892)  72.9 W31-W39 S35-S20
      8 1372  750 (S18W12 L=178)   210     79        0 (1335/ 718-1400/ 779) -59.0 W10-W13 S20-S15
      9 1064  622 (S17E09 L=157)   806    263        0 ( 977/ 508-1144/ 731)  26.6 E17-E01 S21-S12
     10  575  606 (S06E38 L=128)  1061    169        1 ( 526/ 526- 640/ 683) -85.2 E40-E36 S12-N00
     11 1204  458 (S32E06 L=160)   350     68        1 (1178/ 427-1234/ 493)  31.0 E07-E02 S33-S30
     12 1008  346 (S33E23 L=143)  1320    252        2 ( 931/ 265-1097/ 422)  21.3 E32-E14 S37-S29
  

    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