Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1375 1551 (N32W35 L= 33)   469    159        0 (1299/1536-1444/1577)  -9.5 W29-W41 N30-N33
      2 1053 1419 (N30W09 L=  6)   568    130        1 (1002/1373-1121/1457)  16.7 W04-W14 N27-N31
      3 1638 1340 (N14W47 L= 45)  1137    194        3 (1565/1278-1711/1379) -41.6 W41-W52 N09-N17
      4 1019 1145 (N14W01 L=359)   240     71        0 ( 977/1127-1054/1162) -33.7 W00-W03 N12-N15
      5  460 1085 (N17E36 L=321)   777    251        1 ( 359/1038- 540/1159) -29.7 E44-E30 N13-N22
      6 1099  918 (S01W03 L=  0)  1471    362        6 (1049/ 761-1146/1098)  87.3 E01-W04 S11-N10
      7 1404  887 (S08W21 L= 18)   739    247        1 (1312/ 817-1490/ 937) -45.0 W15-W26 S13-S03
      8 1277  829 (S09W12 L= 10)   558    177        0 (1248/ 744-1326/ 909)  52.1 W09-W16 S14-S05
      9 1058  707 (S13E03 L=354)   312    119        0 (1012/ 666-1114/ 745)  29.7 E06-W01 S15-S11
     10 1276  579 (S25W09 L=  6)   507    113        1 (1248/ 528-1305/ 634)  56.3 W06-W11 S28-S22
     11 1112  490 (S29E03 L=354)   228    106        0 (1065/ 476-1173/ 516)  -8.1 E06-W01 S29-S27
     12  935  339 (S37E20 L=337)  1126    290        2 ( 853/ 274- 983/ 439)  40.2 E28-E15 S42-S29
  

    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