Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  680 1686 (N58E18 L=  2)  2043    242        6 ( 580/1617- 789/1728)   0.0 E31-E04 N52-N61
      2 1233 1645 (N41W32 L= 51)  2111    462        2 (1030/1615-1442/1670) -22.8 W15-W47 N34-N47
      3 1622 1522 (N24W56 L= 76)   389     57        0 (1602/1492-1647/1552) -90.0 W54-W57 N21-N25
      4  256 1336 (N36E61 L=319)   662    107        0 ( 219/1286- 297/1367)  35.5 E65-E55 N33-N38
      5  477 1312 (N33E35 L=345)   320     64        1 ( 450/1288- 505/1333)  11.3 E37-E32 N31-N33
      6 1169 1224 (N16W12 L= 32)   862    238        2 (1097/1141-1218/1302) -70.0 W09-W14 N10-N21
      7 1805 1037 (S08W55 L= 74)   370     79        2 (1789/1002-1822/1079)  53.1 W52-W57 S10-S06
      8 1591  689 (S25W32 L= 52)  1368    182        3 (1560/ 615-1636/ 776) -69.0 W30-W36 S31-S18
      9 1288  672 (S19W10 L= 30)   751    191        1 (1245/ 587-1349/ 741)  52.1 W06-W15 S24-S15
  

    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