Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  886 1684 (N50W15 L= 96)   389     85        0 ( 864/1647- 921/1713)  14.0 W11-W19 N47-N51
      2 1654 1645 (N20W70 L=152)  1376    106        4 (1609/1613-1699/1679) -32.0 W65-W77 N17-N23
      3 1673 1575 (N16W65 L=146)   463     72        0 (1648/1545-1698/1609) -90.0 W63-W67 N13-N18
      4 1027 1419 (N28W12 L= 93)   323    117        0 ( 971/1385-1072/1452) -56.3 W09-W13 N25-N31
      5  470 1363 (N41E30 L= 52)   611    165        0 ( 447/1297- 501/1443) -70.0 E32-E27 N35-N46
      6  784 1258 (N26E08 L= 74)   634    208        2 ( 681/1218- 856/1285) -35.0 E14-E03 N22-N30
      7 1870 1150 (S12W63 L=144)  1850    181        1 (1798/1088-1946/1191)  -4.1 W52-W80 S15-S10
      8  444 1135 (N27E34 L= 47)   452     94        3 ( 398/1109- 495/1166)  -7.1 E38-E30 N25-N28
      9  535 1076 (N21E29 L= 53)   338     75        3 ( 514/1037- 562/1107)  63.4 E30-E26 N19-N23
     10  298  986 (N21E50 L= 32)   557    159        1 ( 235/ 941- 350/1016) -31.0 E56-E45 N17-N23
     11  485  874 (N11E37 L= 45)   635    101        2 ( 460/ 828- 520/ 924)  80.5 E38-E33 N08-N14
     12 1523  759 (S23W23 L=104)  1003    391        0 (1380/ 652-1637/ 898)  15.3 W11-W33 S25-S18
     13 1362  798 (S15W13 L= 95)   301     50        1 (1336/ 772-1392/ 818)   0.0 W11-W15 S16-S14
     14  459  793 (N07E41 L= 41)   242     47        1 ( 434/ 768- 485/ 814)  14.0 E42-E38 N05-N07
     15 1785  755 (S33W46 L=127)   704    120        1 (1740/ 713-1821/ 803)  11.3 W40-W50 S34-S30
     16 1217  619 (S22E00 L= 81)   429    203        0 (1159/ 554-1280/ 700)  31.0 E05-W05 S24-S18
     17 1474  622 (S30W16 L= 98)   238     52        0 (1442/ 601-1500/ 639)   0.0 W13-W18 S30-S28
     18  745  426 (S21E35 L= 46)   636    104        2 ( 702/ 389- 794/ 457)   7.1 E39-E31 S22-S19
     19  658  276 (S30E52 L= 29)   734    145        0 ( 611/ 249- 693/ 334)  21.0 E59-E46 S31-S25
  

    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