Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  556 1787 (N31E67 L= 75)  1677    115        1 ( 530/1724- 585/1840) -29.1 E79-E59 N26-N36
      2 1494 1705 (N56W17 L=159)  1569    132        3 (1440/1661-1539/1754)  40.2 W11-W24 N49-N61
      3 1198 1660 (N40E08 L=134)   567    128        1 (1138/1628-1248/1680)   0.0 E12-E03 N37-N41
      4  365 1547 (N11E62 L= 80)  1580    216        2 ( 302/1487- 418/1617)  73.3 E67-E58 N06-N16
      5  891 1483 (N19E20 L=122)  1225    262        1 ( 835/1391- 955/1600) -33.7 E28-E13 N14-N25
      6 1606 1573 (N48W30 L=172)   391     36        0 (1576/1557-1627/1590)  33.7 W26-W32 N46-N50
      7  401 1494 (N09E55 L= 87)  1098    157        2 ( 355/1451- 445/1563)  76.0 E58-E52 N05-N13
      8 1324 1531 (N35W05 L=147)   375     58        1 (1297/1509-1358/1554) -11.3 W02-W07 N32-N35
      9 1653 1169 (N21W36 L=178)   456    118        0 (1599/1149-1703/1190)   6.3 W31-W40 N20-N22
     10 1646  984 (N10W39 L=181)   263     48        0 (1631/ 956-1662/1010)  90.0 W38-W40 N08-N12
     11  391  835 (S30E38 L=104)   260     62        0 ( 363/ 819- 420/ 856)  -9.5 E41-E35 S31-S29
     12  706  812 (S24E14 L=128)   315    113        0 ( 663/ 797- 764/ 834)  26.6 E16-E10 S25-S21
     13  733  764 (S26E11 L=131)   646    328        2 ( 630/ 723- 864/ 800)  29.7 E17-E02 S31-S21
     14 1415  575 (S19W37 L=179)   426     61        1 (1394/ 549-1439/ 606) -26.6 W35-W39 S20-S16
     15 1627  581 (S12W53 L=195)   399     36        1 (1598/ 560-1652/ 603)  33.7 W50-W55 S13-S10
     16  864  535 (S37W04 L=146)   334     59        0 ( 833/ 517- 898/ 558)  45.0 W02-W05 S38-S34
     17 1477  429 (S25W51 L=192)   516     81        0 (1433/ 409-1517/ 449)   6.3 W45-W54 S26-S23
     18  950  293 (S52W23 L=165)   515     85        1 ( 914/ 266- 974/ 328)  68.2 W20-W26 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