Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1596 1515 (N15W52 L=192)   321     74        1 (1559/1495-1630/1530) -31.0 W50-W55 N14-N17
      2 1733 1431 (N07W59 L=199)   583    106        2 (1686/1404-1760/1470)   0.0 W53-W64 N05-N08
      3 1078 1402 (N25W14 L=154)   216     58        1 (1041/1383-1106/1423) -53.1 W12-W15 N23-N27
      4 1869 1190 (S11W62 L=202)  3438    339        9 (1806/1054-1922/1292)  61.7 W55-W73 S19-S05
      5 1013 1252 (N19W06 L=145)   428    127        0 ( 943/1236-1075/1268) -26.6 W01-W09 N16-N21
      6  848 1158 (N18E06 L=133)   312    176        1 ( 779/1106- 914/1201) -53.1 E09-E03 N13-N21
      7 1692 1156 (S06W44 L=183)   321     72        3 (1673/1127-1718/1185) -76.0 W42-W45 S08-S04
      8 1826 1043 (S17W53 L=192)  2345    316        7 (1779/ 868-1866/1155) -87.0 W49-W57 S28-S09
      9  278  975 (N21E51 L= 89)  2056    391        2 ( 180/ 833- 362/1151) -54.5 E62-E45 N12-N33
     10 1674 1044 (S11W39 L=179)   365    140        0 (1634/1007-1724/1085)   0.0 W35-W43 S12-S09
     11  500  974 (N17E33 L=107)   293    114        0 ( 445/ 944- 550/1011) -45.0 E35-E29 N13-N19
     12 1227  804 (S11W05 L=145)   648    342        1 (1179/ 646-1247/ 947)  56.3 E01-W09 S18-S03
     13 1581  889 (S17W29 L=168)   294     97        0 (1532/ 876-1626/ 900) -33.7 W25-W31 S19-S14
     14 1369  562 (S30W07 L=147)  2991    581        8 (1217/ 367-1530/ 782)  31.8 E08-W22 S38-S20
     15 1110  677 (S15E05 L=135)   193    107        0 (1063/ 659-1159/ 691) -38.7 E07-E02 S17-S13
     16  386  616 (S01E52 L= 88)  1275    194        4 ( 305/ 571- 463/ 680) -18.4 E60-E46 S05-N02
     17  559  611 (S05E39 L=101)   880    217        5 ( 505/ 499- 594/ 675)  65.6 E44-E34 S11-N00
     18  757  451 (S19E33 L=107)  1021    328        2 ( 735/ 309- 793/ 594)  46.7 E41-E25 S29-S11
     19  948  501 (S21E20 L=120)   296    117        1 ( 930/ 449- 960/ 570)  54.5 E22-E16 S24-S17
  

    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