Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  475 1665 (N34E55 L=  1)  2025    503        6 ( 432/1592- 542/1744) -25.3 E65-E45 N30-N39
      2  907 1582 (N35E17 L= 39)   374     89        2 ( 863/1557- 948/1611) -15.9 E20-E13 N33-N36
      3 1396 1509 (N36W21 L= 77)   967    175        5 (1343/1384-1441/1604) -68.2 W17-W24 N28-N43
      4  712 1375 (N18E25 L= 31)  1779    502        7 ( 583/1215- 823/1526) -36.5 E38-E14 N09-N26
      5 1090 1436 (N27E01 L= 55)   910    216        2 (1018/1339-1177/1514) -36.9 E07-W05 N21-N31
      6  205 1303 (N07E67 L=348)   660    128        1 ( 186/1270- 235/1336) -76.0 E71-E63 N05-N09
      7  243 1235 (N04E60 L=356)   398    105        0 ( 223/1193- 255/1290)  80.5 E61-E58 N01-N07
      8 1688 1120 (N14W44 L=100)  1265    333        5 (1656/ 981-1723/1285) -81.0 W41-W47 N05-N24
      9 1525 1217 (N18W31 L= 86)   455    167        2 (1482/1149-1553/1264) -66.8 W27-W32 N14-N21
     10 1309 1143 (N11W16 L= 72)   305    112        0 (1258/1105-1357/1186) -29.7 W12-W19 N09-N13
     11 1881 1071 (N13W65 L=120)  2118    230        3 (1852/ 965-1915/1157) -76.0 W60-W70 N06-N18
     12 1093 1068 (N04W04 L= 59)   605    185        1 ( 993/1047-1177/1095)  21.8 E02-W08 N01-N06
     13  374  943 (S12E43 L= 13)   884    213        1 ( 298/ 904- 451/ 995) -13.0 E49-E36 S13-S09
     14  568  942 (S10E28 L= 28)   402    141        1 ( 498/ 923- 627/ 958)  20.6 E32-E23 S11-S08
     15 1180  848 (S08W12 L= 67)   251    125        0 (1144/ 780-1221/ 909)  71.6 W10-W13 S12-S03
     16  253  775 (S24E54 L=  1)   527    131        4 ( 220/ 752- 275/ 809)  39.8 E58-E51 S26-S21
     17 1728  704 (S11W55 L=110)   759    183        1 (1696/ 669-1777/ 758) -35.5 W50-W60 S12-S07
     18 1457  374 (S36W44 L=100)   492     89        1 (1425/ 357-1493/ 400)  21.8 W41-W48 S37-S34
  

    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