Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1137 1871 (N56W40 L=349)   683     60        0 (1122/1826-1158/1908)  28.1 W33-W48 N51-N59
      2 1088 1647 (N38W18 L=327)   918    162        1 (1006/1615-1168/1668)  -9.5 W12-W24 N36-N40
      3 1428 1578 (N27W39 L=348)   753    241        0 (1344/1504-1527/1639) -59.7 W35-W43 N20-N32
      4 1350 1499 (N23W30 L=340)   463    111        1 (1305/1473-1405/1524)  -8.1 W27-W34 N21-N25
      5  781 1412 (N29E09 L=300)   941    383        0 ( 608/1391- 969/1430) -16.9 E20-W03 N24-N31
      6 1188 1339 (N17W16 L=325)   512    226        0 (1090/1319-1296/1355) -13.0 W09-W22 N14-N19
      7  259 1288 (N30E53 L=256)   584    114        1 ( 227/1254- 303/1339) -37.9 E58-E47 N27-N34
      8  770 1293 (N21E11 L=298)   222    127        0 ( 701/1276- 825/1307) -20.6 E15-E07 N20-N23
      9  752 1177 (N15E14 L=295)   821    373        2 ( 616/1080- 887/1259) -43.0 E21-E06 N06-N21
     10 1852 1024 (S14W59 L=  8)   591    108        2 (1830/ 982-1879/1076) -63.4 W56-W61 S16-S10
     11  548  893 (N02E31 L=278)   342    122        0 ( 497/ 852- 598/ 935)  20.6 E35-E27 N00-N03
     12  837  859 (S05E14 L=295)   216    100        0 ( 829/ 804- 850/ 903)  63.4 E15-E12 S08-S02
     13 1901  809 (S28W68 L= 17)  1059     91        0 (1862/ 774-1929/ 855)  24.0 W61-W73 S29-S25
     14 1151  741 (S17W03 L=312)   473    103        1 (1107/ 703-1186/ 784)  31.0 W00-W05 S19-S15
     15  757  682 (S14E22 L=287)   260    113        1 ( 717/ 641- 796/ 722) -63.4 E24-E20 S17-S11
     16  874  442 (S31E22 L=287)   900    251        1 ( 847/ 333- 918/ 558)  49.4 E28-E15 S39-S24
     17 1391  508 (S37W16 L=325)   282    121        0 (1332/ 479-1435/ 539)  12.5 W10-W19 S38-S35
  

    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