Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1022 1726 (N36E18 L= 88)   817    143        2 ( 952/1686-1087/1762)  -9.5 E24-E12 N34-N37
      2  387 1541 (N12E55 L= 51)   686     93        5 ( 359/1498- 427/1582)  80.5 E58-E52 N09-N15
      3 1227 1509 (N26W02 L=108)   361    112        1 (1186/1467-1270/1559)  76.0 W00-W03 N22-N30
      4 1493 1514 (N33W20 L=126)  1406    392        2 (1290/1464-1637/1549)  10.1 W04-W32 N29-N37
      5 1069 1260 (N07E02 L=104)   276     90        0 (1027/1242-1119/1277)  36.9 E04-W00 N05-N08
      6  195 1217 (S09E63 L= 44)   490    101        1 ( 173/1165- 220/1260)  63.4 E64-E60 S12-S06
      7  294 1142 (S13E51 L= 56)   355     92        0 ( 257/1122- 328/1156)  18.4 E53-E47 S14-S11
      8 1825 1095 (N16W52 L=158)   776    112        0 (1781/1046-1861/1137) -18.4 W46-W56 N14-N17
      9 1432 1046 (N03W23 L=129)   309    101        0 (1421/ 993-1454/1092) -59.0 W21-W25 N00-N05
     10  933  990 (S11E05 L=102)   721    355        1 ( 791/ 931-1074/1053)  31.0 E12-W03 S15-S05
     11 1837  837 (N02W59 L=165)  1711    350        3 (1775/ 681-1868/1003) -73.6 W55-W62 S06-N11
     12  300  955 (S24E48 L= 58)   395     65        2 ( 265/ 928- 331/ 981)  18.4 E51-E44 S25-S21
     13 1129  841 (S15W10 L=116)   187     46        0 (1109/ 812-1151/ 865)  76.0 W08-W10 S17-S13
     14 1291  730 (S18W22 L=129)   229     36        1 (1272/ 708-1316/ 753)  71.6 W20-W23 S19-S16
     15 1162  603 (S28W19 L=125)   197     34        0 (1138/ 585-1188/ 619)  56.3 W17-W20 S29-S26
     16 1743  560 (S15W63 L=170)   552    116        2 (1719/ 505-1766/ 605) -56.3 W61-W66 S17-S11
  

    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