Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1112 1901 (N59E37 L=195)   586     40        1 (1087/1880-1131/1924) -45.0 E43-E33 N55-N63
      2 1198 1844 (N56E21 L=211)  1424    145        4 (1146/1797-1265/1890) -14.7 E31-E12 N52-N60
      3 1415 1754 (N58W05 L=237)   400     47        1 (1395/1729-1442/1771)  60.3 W03-W08 N53-N60
      4  883 1627 (N27E27 L=205)  1079    303        1 ( 762/1539- 980/1693) -10.8 E38-E17 N24-N28
      5 1300 1503 (N32W03 L=236)   781    146        2 (1247/1446-1348/1560) -26.6 E01-W07 N29-N34
      6  970 1409 (N16E14 L=219)   304    124        0 ( 924/1369-1022/1433)   0.0 E16-E09 N14-N17
      7  229 1274 (S09E63 L=169)  1169    120        2 ( 195/1220- 270/1316)  49.4 E67-E59 S12-S05
      8 1272 1233 (N14W09 L=241)   675    178        0 (1249/1148-1305/1309) -79.7 W07-W10 N08-N19
      9  592  988 (S17E25 L=207)   406    113        1 ( 530/ 972- 642/1002)  14.0 E29-E21 S18-S15
     10 1010  776 (S19W06 L=238)   377    119        2 ( 954/ 757-1063/ 806)   8.1 W02-W09 S19-S17
     11  909  703 (S25W02 L=234)   241    145        0 ( 853/ 659- 951/ 757)  76.0 W00-W03 S29-S21
     12 1446  649 (S13W36 L=269)   344    111        1 (1410/ 599-1468/ 692) -33.7 W32-W39 S15-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