Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  939 1811 (N40E37 L=347)   621     65        0 ( 901/1787- 972/1836)  31.0 E39-E32 N37-N42
      2  540 1592 (N14E49 L=335)  2166    337        0 ( 494/1458- 581/1754) -56.3 E60-E43 N04-N25
      3  462 1632 (N15E58 L=326)  1123    186        1 ( 436/1548- 504/1714) -56.3 E63-E53 N08-N20
      4 1353 1568 (N36W03 L= 27)   436     60        1 (1334/1538-1372/1607) -45.0 W01-W05 N33-N38
      5  607 1475 (N08E38 L=346)  1949    368        2 ( 475/1397- 745/1568)  55.5 E45-E32 N00-N16
      6 1373 1489 (N31W07 L= 31)   412     66        2 (1340/1463-1406/1519)  11.3 W04-W09 N29-N33
      7 1030 1399 (N15E10 L= 14)   301    124        0 ( 990/1349-1074/1443) -23.2 E13-E06 N13-N16
      8  207 1278 (S11E66 L=318)  1579    172        4 ( 124/1246- 258/1326)  12.8 E83-E60 S13-S07
      9  192 1170 (S18E64 L=320)   490     92        2 ( 172/1131- 208/1223) -80.5 E66-E62 S20-S14
     10 1316 1150 (N09W13 L= 36)   255     56        1 (1301/1116-1329/1183) -56.3 W11-W13 N07-N10
     11  418 1130 (S16E40 L=343)   305     63        0 ( 392/1107- 440/1162)  63.4 E42-E38 S17-S13
     12  174  960 (S31E64 L=320)  1732    188        1 ( 118/ 888- 230/1021)  23.2 E77-E56 S35-S26
     13  726  786 (S28E11 L= 13)   881    239        1 ( 618/ 747- 834/ 832)  39.8 E17-E05 S32-S22
     14 1378  717 (S12W30 L= 53)   292    148        1 (1310/ 692-1449/ 752)  12.5 W24-W34 S13-S11
     15  294  589 (S54E54 L=330)   419     59        0 ( 270/ 573- 327/ 610)   0.0 E59-E47 S54-S52
     16 1219  580 (S25W26 L= 50)   323    122        1 (1156/ 570-1283/ 591)  29.7 W21-W29 S26-S22
     17 1495  530 (S18W46 L= 70)   990    170        2 (1428/ 505-1581/ 571)  48.8 W40-W50 S21-S13
     18 1097  409 (S39W28 L= 52)  1849    359        3 ( 958/ 326-1239/ 485)  68.2 W22-W32 S49-S29
     19  928  382 (S47W19 L= 42)   291     27        1 ( 908/ 360- 947/ 407) -63.4 W16-W20 S48-S44
  

    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