Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  521 1424 (N03E41 L=130)   282     61        2 ( 506/1395- 538/1453)  90.0 E43-E40 N01-N05
      2  815 1353 (N06E21 L=151)   339    134        1 ( 752/1316- 885/1383)  45.0 E23-E17 N02-N09
      3 1338 1321 (N18W10 L=182)   509    206        2 (1259/1237-1415/1364) -14.0 W04-W16 N15-N19
      4  291 1229 (S12E53 L=119)   477    189        1 ( 241/1182- 331/1280) -18.4 E60-E48 S14-S10
      5 1260 1120 (N05W10 L=182)   359    130        1 (1198/1103-1332/1138)  26.6 W06-W14 N03-N07
      6  455  886 (S29E32 L=140)   806    239        2 ( 369/ 835- 544/ 949)  42.5 E38-E26 S34-S23
      7 1634  705 (S06W47 L=218)   624    158        1 (1582/ 682-1696/ 728)  11.3 W42-W52 S07-S04
      8 1689  578 (S10W57 L=229)  1039    276        5 (1616/ 544-1767/ 625)  42.0 W52-W65 S14-S05
      9 1008  512 (S36W16 L=188)  1196    261        4 ( 893/ 440-1106/ 585)   0.0 W06-W25 S38-S33
     10 1335  549 (S24W34 L=206)   296     89        0 (1284/ 530-1380/ 573)  51.3 W31-W35 S26-S20
     11 1278  454 (S31W35 L=207)   315    100        0 (1250/ 421-1309/ 498) -23.2 W31-W38 S31-S28
  

    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