Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  957 1469 (N33E03 L= 29)   340     61        0 ( 934/1433- 987/1500) -53.1 E04-W00 N30-N34
      2  252 1431 (N30E67 L=325)  2895    136        7 ( 211/1358- 289/1487) -22.8 E79-E58 N25-N33
      3  749 1300 (N22E17 L= 15)   256     79        1 ( 724/1277- 779/1323) -33.7 E18-E15 N20-N23
      4  483 1270 (N20E37 L=355)  1057    228        1 ( 390/1241- 586/1303)   0.0 E45-E28 N18-N22
      5 1071 1250 (N18W04 L= 36)   202     86        0 (1031/1238-1107/1275) -21.8 W01-W06 N17-N19
      6  835 1162 (N13E11 L= 20)   301     62        0 ( 795/1147- 871/1184)   0.0 E13-E09 N12-N14
      7 1597  970 (S02W37 L= 69)   787    235        0 (1576/ 864-1620/1082) -90.0 W35-W38 S08-N04
      8  596  971 (N02E27 L=  5)   318    131        0 ( 542/ 946- 659/ 994) -15.9 E30-E23 N00-N03
      9 1224  954 (S01W12 L= 44)   287     63        0 (1208/ 917-1240/ 993)  76.0 W10-W13 S03-N01
     10  783  817 (S08E16 L= 16)   365    130        0 ( 731/ 778- 838/ 842)  29.7 E19-E12 S10-S06
     11 1513  799 (S12W31 L= 63)   700    208        0 (1428/ 747-1597/ 836)  18.4 W25-W37 S14-S10
  

    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