Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  655 1680 (N46E36 L=128)  2876    373        4 ( 533/1600- 792/1734)  10.3 E45-E22 N38-N50
      2 1285 1632 (N44W20 L=184)   396     69        0 (1269/1592-1305/1662)  59.0 W18-W22 N41-N46
      3  992 1614 (N42E05 L=159)   384     66        2 ( 964/1587-1018/1643)  68.2 E07-E02 N39-N44
      4 1411 1540 (N37W29 L=193)   281     55        0 (1384/1526-1442/1561)  11.3 W26-W31 N36-N38
      5  962 1409 (N27E05 L=159)   215     81        0 ( 916/1395- 993/1431)  11.3 E08-E03 N26-N28
      6  518 1231 (N14E35 L=130)   519    126        2 ( 466/1189- 549/1264) -26.6 E38-E32 N11-N15
      7  599 1247 (N15E29 L=135)   265     91        0 ( 559/1232- 645/1258)   9.5 E31-E25 N14-N15
      8  637 1191 (N12E26 L=139)   341     61        1 ( 620/1153- 665/1227) -78.7 E26-E23 N09-N14
      9  516 1152 (N09E34 L=130)   786    259        2 ( 447/1107- 627/1215) -21.0 E39-E25 N06-N12
     10 1600  575 (S25W44 L=208)   499    121        0 (1560/ 546-1637/ 603) -20.6 W39-W47 S26-S22
     11 1422  418 (S37W34 L=198)  1243    220        2 (1358/ 386-1510/ 468)  -4.4 W27-W40 S39-S32
     12 1008  389 (S40W01 L=165)  1198    300        2 ( 874/ 350-1161/ 421)  15.3 E10-W12 S43-S36
  

    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