Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  833 1648 (N46E17 L=342)   517    106        2 ( 785/1615- 881/1680) -15.3 E22-E11 N42-N48
      2  956 1591 (N41E05 L=353)   257     22        0 ( 934/1575- 980/1604)  18.4 E07-E03 N39-N42
      3 1519 1391 (N26W36 L= 35)   528    121        1 (1498/1335-1544/1449)  81.9 W34-W38 N22-N29
      4  616 1271 (N19E27 L=331)   924    288        0 ( 552/1163- 698/1382) -45.0 E33-E20 N12-N25
      5 1248 1159 (N12W14 L= 12)   946    273        0 (1169/1087-1363/1237) -36.9 W09-W21 N07-N16
      6 1015 1117 (N09E01 L=358)   219    125        0 ( 992/1048-1033/1166) -81.9 E01-W00 N05-N12
      7 1717 1042 (N03W48 L= 46)   354     64        0 (1696/1008-1733/1079)  53.1 W46-W49 N01-N05
      8 1567  895 (S05W36 L= 34)   259     74        0 (1553/ 858-1583/ 928) -76.0 W34-W37 S07-S03
      9 1107  793 (S11W05 L=  3)   546    184        1 (1053/ 755-1185/ 847)  20.6 W01-W09 S13-S07
     10  996  738 (S14E02 L=356)   208     96        0 ( 973/ 688-1014/ 776)  71.6 E03-W00 S17-S11
     11  901  670 (S19E08 L=350)   245     54        1 ( 882/ 642- 922/ 698)  53.1 E09-E06 S20-S16
     12 1632  596 (S25W46 L= 44)   317     96        1 (1599/ 571-1679/ 628)  31.0 W42-W49 S26-S23
  

    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