Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1574 1686 (N26W66 L= 60)   544    129        1 (1515/1663-1637/1698) -14.0 W59-W75 N23-N27
      2 1405 1576 (N26W43 L= 37)   486    109        3 (1375/1519-1433/1622)  66.8 W40-W45 N22-N29
      3 1360 1405 (N18W31 L= 25)   384    143        0 (1326/1334-1402/1452)  32.0 W26-W35 N14-N20
      4 1139 1417 (N25W19 L= 12)   346    111        0 (1084/1394-1196/1446) -50.2 W15-W21 N22-N28
      5 1018 1313 (N23W08 L=  2)   526    190        0 ( 939/1268-1090/1344) -45.0 W03-W10 N18-N26
      6  737 1207 (N24E12 L=342)   314    112        0 ( 709/1154- 762/1252) -81.9 E13-E11 N20-N27
      7  287 1208 (N35E50 L=304)   744    141        0 ( 237/1169- 355/1250) -23.2 E56-E42 N32-N38
      8  763 1126 (N19E13 L=341)   213    117        0 ( 719/1095- 814/1161) -56.3 E14-E10 N16-N22
      9 1907 1093 (S17W67 L= 60)  1146    139        1 (1851/1052-1957/1117) -21.8 W58-W78 S21-S13
     10 1654 1035 (S11W38 L= 32)   632    173        0 (1616/ 958-1707/1096)  26.6 W33-W44 S13-S08
     11 1851  987 (S20W56 L= 50)   470     79        2 (1836/ 951-1868/1026) -80.5 W54-W57 S23-S17
     12 1271  782 (S13W08 L=  1)   868    300        1 (1211/ 667-1326/ 918)  40.2 W01-W14 S18-S07
     13  516  844 (N09E35 L=319)   725    171        2 ( 465/ 779- 580/ 916) -66.0 E37-E32 N04-N13
     14 1588  823 (S21W29 L= 23)   342     68        1 (1571/ 791-1604/ 861)  76.0 W27-W29 S22-S18
     15 1290  531 (S30W02 L=356)   414     40        1 (1268/ 504-1317/ 562)  14.0 W00-W04 S31-S27
     16 1031  535 (S22E13 L=341)   315    112        0 ( 987/ 518-1089/ 559) -26.6 E15-E09 S23-S19
     17  838  433 (S23E29 L=325)   498    182        0 ( 794/ 361- 868/ 495)  35.0 E35-E25 S26-S19
     18  820  230 (S37E44 L=309)   448     84        0 ( 797/ 194- 844/ 274)  49.4 E47-E41 S40-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