Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1045 1430 (N27W13 L=312)   253     68        0 (1015/1404-1083/1453)   0.0 W10-W16 N25-N27
      2 1677 1406 (N07W52 L=351)   483    127        1 (1624/1380-1720/1434) -18.4 W48-W57 N05-N09
      3  180 1307 (N43E65 L=234)   614     93        1 ( 137/1290- 226/1322)  -3.2 E74-E56 N41-N43
      4 1136 1158 (N10W10 L=308)  1325    531        1 ( 930/1031-1331/1245) -46.6 W00-W17 S02-N17
      5 1584 1185 (S01W37 L=335)   764    204        0 (1516/1135-1640/1220) -16.7 W30-W40 S04-N00
      6 1649 1082 (S09W38 L=337)   243     37        1 (1630/1049-1661/1107)  71.6 W37-W39 S10-S07
      7  608 1012 (N15E24 L=275)   820    183        1 ( 573/ 943- 662/1102) -84.3 E25-E22 N10-N20
      8 1603  774 (S25W28 L=327)  1885    415        1 (1509/ 614-1734/ 931)  28.6 W18-W41 S32-S19
      9  515  807 (N06E35 L=264)   665    115        2 ( 462/ 784- 581/ 842) -29.7 E38-E31 N03-N08
     10 1168  615 (S21E03 L=296)   627    167        2 (1117/ 553-1212/ 685)  32.0 E07-W01 S23-S18
     11  545  528 (S11E43 L=256)  1634    368        2 ( 476/ 465- 618/ 616) -79.7 E47-E37 S15-S04
     12 1431  498 (S37W10 L=309)   365    108        0 (1389/ 469-1467/ 541)  15.9 W07-W14 S38-S35
     13  722  452 (S19E35 L=264)   293    121        0 ( 696/ 409- 755/ 497)  63.4 E36-E32 S22-S16
     14  503  393 (S18E54 L=244)  1691    302        4 ( 384/ 329- 601/ 496) -60.3 E60-E47 S23-S09
     15 1019  451 (S27E17 L=282)   648     94        2 ( 978/ 414-1069/ 492)   8.1 E20-E13 S29-S24
     16 1289  326 (S45E05 L=293)   278     58        0 (1272/ 298-1303/ 360)  38.7 E07-E02 S47-S42
     17  852  254 (S36E39 L=260)   786     96        1 ( 815/ 224- 892/ 284)  12.5 E43-E34 S38-S33
  

    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