Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  537 1613 (N35E46 L= 66)   603    195        1 ( 456/1559- 598/1669)  55.0 E51-E43 N30-N40
      2 1333 1540 (N37W19 L=131)   626    185        0 (1253/1483-1393/1592) -32.5 W12-W23 N33-N41
      3  802 1466 (N28E19 L= 93)  1040    203        4 ( 741/1385- 870/1542) -42.3 E25-E13 N22-N32
      4 1096 1368 (N24W02 L=114)   382    183        1 (1069/1292-1117/1447) -84.3 W00-W03 N18-N28
      5 1399 1338 (N24W23 L=135)   306    110        1 (1379/1290-1419/1399) -90.0 W21-W24 N21-N28
      6 1818 1158 (N15W58 L=170)   428    106        2 (1800/1105-1831/1199) -90.0 W56-W60 N11-N17
      7  249 1069 (S02E56 L= 56)   486     88        1 ( 234/1017- 263/1107) -68.2 E57-E54 S05-N00
      8 1710 1021 (N06W47 L=159)   453    101        1 (1683/ 984-1749/1059) -33.7 W44-W50 N04-N08
      9  282  963 (S08E52 L= 60)   347    144        0 ( 236/ 936- 328/ 992)   6.3 E56-E47 S10-S06
     10 1009  701 (S18W02 L=114)  2694    743       10 ( 735/ 531-1240/ 950)  48.0 E16-W13 S32-S01
     11 1659  733 (S12W46 L=158)   416     81        3 (1633/ 707-1685/ 763)  90.0 W44-W49 S13-S09
     12  231  660 (S29E65 L= 47)  1104    118        1 ( 193/ 608- 273/ 690)  14.0 E80-E56 S33-S26
     13 1611  582 (S22W47 L=159)   999    221        7 (1530/ 545-1684/ 628)  28.6 W39-W54 S24-S18
     14 1622  509 (S27W52 L=164)   702    105        3 (1592/ 455-1650/ 544) -56.3 W48-W56 S30-S24
     15  461  451 (S42E44 L= 68)   998    122        3 ( 409/ 407- 501/ 492)  25.2 E55-E38 S46-S38
     16 1691  442 (S31W67 L=179)   883     67        3 (1665/ 419-1717/ 471) -21.8 W61-W73 S33-S29
     17 1608  367 (S38W63 L=175)   911     71        1 (1581/ 329-1633/ 400) -14.7 W56-W75 S40-S34
     18  676  318 (S51E26 L= 86)   426     85        1 ( 648/ 298- 717/ 339)  29.1 E29-E20 S53-S48
  

    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