Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1481 1727 (N48W55 L=184)   448     58        0 (1457/1702-1503/1752)  16.7 W50-W60 N46-N49
      2  877 1420 (N30E08 L=121)   648    183        1 ( 813/1350- 912/1498) -59.0 E13-E05 N25-N35
      3 1279 1336 (N22W19 L=148)   736    243        0 (1203/1243-1344/1371) -45.0 W13-W22 N16-N24
      4  243 1293 (N23E61 L= 68)  1083    122        2 ( 217/1236- 272/1352) -60.3 E63-E57 N19-N26
      5 1816 1270 (N13W63 L=191)   544     75        3 (1785/1236-1839/1302) -76.0 W59-W65 N11-N15
      6  388 1284 (N23E45 L= 84)   283     62        1 ( 356/1265- 419/1297)   0.0 E48-E42 N21-N23
      7  176 1239 (N19E70 L= 59)  1112    108        1 ( 153/1187- 196/1286)  66.8 E73-E66 N15-N22
      8 1554 1197 (N12W36 L=165)   341     76        2 (1530/1167-1587/1228) -45.0 W34-W38 N09-N13
      9 1835  939 (S07W60 L=188)   590     70        2 (1821/ 904-1856/ 978) -68.2 W57-W62 S09-S04
     10 1634  839 (S11W40 L=169)   505     93        1 (1611/ 788-1658/ 889)  63.4 W38-W42 S14-S07
     11 1569  837 (S10W35 L=164)   282     63        4 (1534/ 816-1602/ 858)   0.0 W32-W37 S11-S09
     12 1538  754 (S16W33 L=162)   802    232        5 (1423/ 709-1602/ 816)  -8.7 W24-W37 S18-S11
     13  950  741 (S13E06 L=123)   470    201        0 ( 904/ 661-1016/ 811)  48.4 E09-E01 S17-S08
     14  485  734 (S12E38 L= 91)   304     92        3 ( 462/ 684- 508/ 766)  51.3 E39-E35 S14-S09
     15 1797  691 (S23W59 L=188)  1137    130        1 (1748/ 657-1843/ 736)  16.7 W53-W65 S24-S20
     16 1539  670 (S21W34 L=162)   525    136        2 (1480/ 647-1585/ 697)  -7.1 W29-W37 S22-S19
     17 1512  602 (S26W32 L=161)   365     89        2 (1483/ 571-1535/ 633) -45.0 W29-W34 S27-S23
     18  595  533 (S25E33 L= 95)  1097    145        3 ( 570/ 459- 641/ 600)  59.0 E37-E30 S30-S20
     19  484  433 (S33E48 L= 81)  1001    175        2 ( 408/ 414- 569/ 454)  -3.6 E56-E39 S34-S31
  

    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