Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  789 1708 (N57W05 L= 42)   422     50        1 ( 764/1677- 813/1739)  45.0 W02-W09 N54-N59
      2  447 1545 (N54E37 L=359)   387     47        1 ( 417/1525- 471/1565)   0.0 E41-E34 N52-N56
      3 1206 1356 (N21W20 L= 56)   262    112        0 (1159/1320-1240/1371) -59.0 W17-W21 N18-N23
      4  962 1164 (N17E00 L= 36)   797    190        1 ( 881/1104-1028/1230) -58.0 E03-W03 N12-N22
      5 1784 1168 (S06W54 L= 91)   354     56        2 (1760/1141-1805/1197)  33.7 W52-W55 S07-S03
      6  980  927 (N03E05 L= 32)   973    366        1 ( 906/ 779-1046/1081)  51.8 E11-W02 S04-N10
      7  676 1015 (N14E21 L= 16)   257    101        1 ( 661/ 970- 706/1058)  51.3 E22-E18 N11-N16
      8  303  978 (N19E50 L=346)   390     79        0 ( 290/ 932- 315/1020) -90.0 E51-E48 N16-N21
      9  241  981 (N20E57 L=339)   392     74        0 ( 220/ 950- 263/1015)  31.0 E59-E54 N18-N22
     10 1768  906 (S20W47 L= 83)  1424    259        1 (1717/ 783-1823/1003)  34.7 W40-W53 S26-S15
     11 1415  958 (S06W21 L= 57)   212     44        0 (1399/ 929-1432/ 988)  76.0 W20-W21 S08-S04
     12 1828  836 (S27W53 L= 89)  1649    310        0 (1734/ 745-1911/ 945)  11.8 W41-W66 S31-S22
     13 1135  612 (S20E04 L= 33)   442    178        0 (1059/ 575-1200/ 631) -32.0 E08-W00 S23-S17
     14 1091  451 (S29E11 L= 25)   680    167        3 (1028/ 410-1161/ 505)  15.3 E16-E05 S30-S26
     15 1359  440 (S38W07 L= 43)   918    143        2 (1309/ 393-1404/ 503)  33.7 W02-W11 S41-S34
     16 1146  297 (S43E14 L= 22)   353     39        0 (1125/ 268-1165/ 326)  80.5 E15-E12 S45-S39
  

    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